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Abstract. We report the hrst test of isotropy of the Universe in the matter dom¬ 
inated epoch using the Lyman-a forest data from the high redshift quasars {z > 2) 
from SDSS-III BOSS-DR9 datasets. Using some specihed data cuts, we obtain the 
probability distribution function (PDF) of the Lyman-a forest transmitted flux and 
use the statistical moments of the PDF to address the isotropy of the Universe. In an 
isotropic Universe one would expect the transmitted flux to have consistent statistical 
characteristics in different parts of the sky. We trisect the total survey area of 3275 
deg^ along the galactic latitude and using quadrant convention. We also make three 
subdivisions in the data for three different signal-to-noise-ratios (SNR). Finally we ob¬ 
tain and compare the statistical moments in the mean redshifts of 2.3, 2.6 and 2.9. We 
hnd, that the moments from all patches agree at all redshifts and at all SNRs, within 
3a uncertainties. Since Lyman-a transmitted flux directly maps the neutral hydrogen 
distribution in the inter galactic medium (IGM), our results indicate, within the limited 
survey area and sensitivity of the data, the distribution of the neutral hydrogen in the 
Universe is consistent with isotropic distribution. We should mention that we report 
few deviations from isotropy in the data with low statistical significance. Increase in 
survey area and larger amount of data are needed to make any strong conclusion about 
these deviations. 
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1 Introduction 

Our understanding of the physical cosmology strongly depends on the data that we 
observe and the model assumptions that we make. Assumptions in the standard model 
sometime pose problems in understanding the true nature of the Universe. Statistical 
isotropy is one of such assumptions which we assume in almost all the data analysis. 
In this assumption we claim that the statistical properties of the observables in the 
Universe are same in all direction. While in parameter estimations from the datasets, 
isotropy is assumed, there have been a number of tests to conhrm this assumption, 
using data from different cosmological probes, such as Cosmic Microwave Background 
(CMB) [1-18] and Large Scale Structure (LSS) [19-30]. In this paper we investigate 
isotropy of the matter dominated Universe using the Lyman-a forest datasets in the 
redshift range z ^ 2 —3. Spectra of high redshift quasars contain absorption lines that 
trace the components of the IGM along the line of sight. In the case of Lyman-a forest, 
we hnd absorption lines from the hrst ionization state of the Hydrogen atoms. The 
wavelength of the absorption identihes the redshift of the neutral hydrogen cloud. An¬ 
alyzing these absorption lines, collectively from a number of quasars, we can constrain 
properties of the IGM. In particular, the Lyman-a transmitted flux (F) (absorption 
w.r.t the estimated continuum spectrum) can be related to the dark matter overdensity 
(5) as F = Fexp[-A(l + 5 ) 2 -o-T 7 -i)] [ 31 ]^ where, F is the mean transmitted flux, 7 — 1 
dictates the temperature density relation and A is a redshift dependent constant. We 
shall address the statistics of the transmitted flux in the Lyman-a region at different 
redshifts. To have a theoretical model independent analysis, in this work we only con¬ 
sider the statistical properties of the observational data. In the last decade, detection 
of large number of quasar spectrum with high SNR enables us to perform this type 
of analysis. We use Baryon Oscillation Spectroscopic Survey (BOSS) DR9 from Sloan 
Digital Sky Survey (SDSS)-III [32, 33]. We perform our test in three different redshifts 
in the matter dominated epoch that enables us to track the isotropy along time. In 
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each redshift, the isotropy is tested in medium, high and highest SNR of detection. We 
divide the sky in three patches using two different patch selection criteria. Since SDSS 
is a ground base survey which only has partial sky coverage, our test of isotropy will 
be limited by the survey area. Throughout the analysis we closely follow the survey 
parameters for the selection of data pixels from the quasar spectrum in order to obtain 
the PDF of the transmitted flux. 

The paper is organized as follows : In section 2 we discuss the data from BOSS- 
DR9 and the selection criteria for the data pixels for our analysis. We also outline 
the error estimation procedure for the flux PDF. Afterwards we discuss and tabulate 
the properties of the sky patches selected. In the results section 3 we discuss the main 
outcome of our analysis. Finally in section 4 we conclude along with highlighting future 
prospects. 

2 Data analysis 

We use the latest available BOSS DR9 quasar Lyman-a forest data [34]. The ninth 
data release contains 54,468 spectra of Lyman-a quasars with redshifts more than 2.15. 
They provide Lyman-a forest data in the redshift range z ^ 2 — 5.7. However due to 
a smaller number of quasars with redshift higher than 2 ; ~ 3 we restrict our analysis 
in the redshift z ^ 2 — 3. BOSS-DR9 has a survey area of 3275 square degrees and 
hence our test of isotropy will be limited by this area. As has been discussed in [34], 
out of 87,822 total spectra, the set of 54,468 spectra was selected after removing low 
redshift quasars, quasars having broad absorption lines, too low SNR and negative 
continuum. Note that the continuum to each quasar spectrum is estimated using 
mean-flux regulated principal component analysis (MF-PCA) techniques [35] and they 
are provided along with the data. We would like to mention that even at this stage not 
all the data are appropriate for the isotropy test due to low SNR, damped Lyman-a 
(DBA), limited exposures and few other characteristics. Below we point out the data 
cuts that we have used to select a spectrum. 

2.1 Quasar selection, data cuts 

Our selection of data closely follows the BOSS criteria used to obtain constraints on 
the IGM [36]. As we mentioned earlier, we test the isotropy of the sky in three redshift 
bins with the central redshifts being 2.3, 2.6 and 2.9. The redshift bins have a bin 
width of Az = 0.3. In each of the redshift bins we impose the following data cuts. 
Firstly since large number of low SNR quasars can have a systematic bias in our results, 
we reject all the quasars with SNR < 6. The rest of the quasars are binned in good 
(6 < SNR < 8), better (8 < SNR < 10) and best (SNR > 10) categories. In the quasar 
rest frame we use 1041 — 1185A as the Lyman-a domain. 

Lyman-a forest typically refers to neutral hydrogen atoms with column density 
of 10^^ atoms per cm^ in the line of sight. The presence of DBA in a spectrum in¬ 
dicates extremely high column density (0(10^°) atoms per cm^). Following [36] we 
discard spectra with identihed DBA in the sightlines. The provided data uses a DBA 
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concordance catalogue for the identification of DLA’s. The detection efficiency of 
the catalogue decreases below neutral hydrogen column densities less than 10^°'^/cm^. 
Hence, below this criteria, in our analysis, for each of the forest we correct the flux for 
the damping wings of DLA in the sightlines, provided with the data [34], We leave 
the spectra with neutral hydrogen column density of more than 10^°'^/cm^ out of our 
analysis. 

As mentioned earlier, the survey provides the continuum of each spectrum esti¬ 
mated using MF-PCA technique. Note that the transmitted flux is the observed flux 
w.r.t. the estimated continuum and hence a bad continuum estimation shall bias the 
calculated flux *. Here too, following the survey criteria we reject the quasars which do 
not provide good £t to the spectrum redwards to the Lyman-a forest. We also reject 
spectra which had less than three individual exposures. 

Our hnal selection criteria depends on the resolution of the BOSS spectographs. 
As has been mentioned in [36], BOSS spectographs do not resolve the Lyman-a forest. 
To evade this systematic effect that can reflect in the transmitted flux, we use similar 
procedure used in [36]. We stack the spectra according to the intrinsic wavelength 
dispersion (cxdisp) at the Lyman-a wavelength of each central redshift bin. Note that 
in each pixel the cxdisp is provided along with the spectrum. From each redshift bin we 
then discard the 5% spectra from below and 10% spectra from above. 

Each quasar spectrum contain certain number of data pixels in the Lyman-a 
region. In a redshift bin (say, 2.15 < z < 2.45) some quasars may have all the 
Lyman-a pixels, some will have partial pixels from either blue end or red end of the 
Lyman-a domain. We here consider all the quasars satisfying the above criteria and 
also contributing more than 30 pixels in the redshift window. 

2.2 Selection of sky patches 

Since BOSS-DR9 surveys in 3275 deg^ area, we need to make patches in the sky in 
particular ways. First we convert the coordinate of the quasars from J2000 equatorial 
system to galactic coordinate (/,6) system. Below, in Fig. 1 we provide our selection 
of patches with different colors. In the hrst system we divide the sky using quadrant 
convention. Here we divide the sky in three patches where one patch remains in the 
southern hemisphere (red, b < 0°) and the other two patches are in northern hemisphere 
{b > 0°). The two patches in the northern hemisphere are divided in I > 180° (blue) 
and I < 180° (green) parts. In principle, the quadrant convention divides the southern 
hemisphere into two parts (/ > 180° and I < 180°) too, but since there are barely any 
data in 6 < 0°,/ > 180°, we take only one patch from southern hemisphere. In our 
other selection, we divide the sky according to galactic latitude. The red patch remains 
same as it stays in the southern hemisphere. The patches in the northern hemisphere 
are divided in one patch with b > 50° (magenta) and another with 0° < 6 < 50° 
(black). Our hrst selection of patches is provided in the left of Fig. 1 and the other 
selection is shown to the right. 

In the following table 1 we tabulate the essential information about the selected 
patches and the overall sample in different redshift and different SNR. For different 

*The effects of different continuum estimations are discussed in [d7] 
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Figure 1. The B0SS-DR9 survey area in galactic coordinates and our selection of sky patches for 
the test of isotropy. To accommodate comparable number of quasars in each patch we have divided 
the sky in two manner, namely patch selection 1 (left) and patch selection 2 (right). 


patches we provide the number of quasars and the number of pixels that pass all our 
data cuts discussed in the paragraph before. Note that apart from the red patch 
that stays in the southern hemisphere, other patches contain quasar number and the 
pixels that are comparable to each other. As we go towards higher redshifts we hud 
lesser number of data pixels which is expected as we know there are not many high 
redshift quasars. This decrease in number of high redshift quasars limit our analysis 
to z < 3. However, in each SNR bin we have enough pixels to perform a relatively 
robust statistical analysis. Also note that compared to 6 < SNR < 8 bin, in the two 
higher SNR bins the numbers of quasars or pixels are approximately half. We should 
mention that this is not a concern since we assume the data from different SNR are 
independent and we compare the statistical properties of the transmitted flux in each 
SNR separately. Moreover the lack of large number of quasars in higher SNR bins are 
compensated by the less dispersion in the data from higher SNR quasar samples. 
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Redshift range( 2 ) 

SNR 

Patch location 

Number of quasars 

Number of Lyman-a pixels 



Complete sky 

1091 

285878 



b<0° (Red) 

226 

60631 



b>0°, 1 < 180° (Green) 

390 

100289 


6-8 

fe > 0°, l> 180° (Blue) 

475 

124958 



0° < fe < 50° (Black) 

401 

106772 



b > 50° (Magenta) 

464 

118475 



Complete sky 

498 

130242 



6 < 0° (Red) 

106 

28686 



fe > 0°, 1 < 180° (Green) 

185 

48158 

2.15 - 2.45 (5 = 2.3) 

8-10 

fe > 0°, l> 180° (Blue) 

207 

53398 



0° < fe < 50° (Black) 

168 

45290 



b > 50° (Magenta) 

224 

56266 



Complete sky 

567 

147679 



6 < 0° (Red) 

122 

30899 



fe > 0°, 1 < 180° (Green) 

187 

47806 


> 10 

fe > 0°, l> 180° (Blue) 

258 

68974 



0° <b< 50° (Black) 

228 

62655 



b > 50° (Magenta) 

217 

54125 



Complete sky 

975 

230165 



6 < 0° (Red) 

221 

53062 



6 > 0°, 1 < 180° (Green) 

373 

86995 


6-8 

fe > 0°, l> 180° (Blue) 

381 

90108 



0° < fe < 50° (Black) 

348 

80712 



b > 50° (Magenta) 

406 

96391 



Complete sky 

475 

107539 



6 < 0° (Red) 

no 

25099 



6 > 0°, 1 < 180° (Green) 

181 

41564 

2.45 - 2.75(z = 2.6) 

8-10 

b> 0°, l> 180° (Blue) 

184 

40876 



0° < fe < 50° (Black) 

141 

32220 



b > 50° (Magenta) 

224 

50220 



Complete sky 

607 

144797 



6 < 0° (Red) 

139 

33690 



6 > 0°, 1 < 180° (Green) 

224 

54156 


> 10 

fe > 0°, l> 180° (Blue) 

244 

56951 



0° < fe < 50° (Black) 

248 

56434 



b > 50° (Magenta) 

220 

54673 



Complete sky 

628 

139373 



6 < 0° (Red) 

135 

30291 



fe > 0°, 1 < 180° (Green) 

250 

57698 


6-8 

fe > 0°, l> 180° (Blue) 

243 

51384 



0° < fe < 50° (Black) 

225 

50412 



b > 50° (Magenta) 

268 

58670 



Complete sky 

372 

87257 



6 < 0° (Red) 

75 

17879 



fe > 0°, 1 < 180° (Green) 

150 

35068 

2.75 - 3.05 (2 = 2.9) 

8-10 

fe > 0°, l> 180° (Blue) 

147 

34310 



0° < 6 < 50° (Black) 

112 

25351 



b > 50° (Magenta) 

185 

44027 



Complete sky 

432 

97826 



6 < 0° (Red) 

111 

25645 



fe > 0°, 1 < 180° (Green) 

174 

37713 


> 10 

fe > 0°, l> 180° (Blue) 

147 

34468 



0° < 6 < 50° (Black) 

164 

39660 



b > 50° (Magenta) 

157 

32521 


Table 1. The number of quasars and the number of Lyman-a pixels that contribute in 
different bins and patches. We bin our samples in three different redshift bins with mean 
redshift being 2.3, 2.6 and 2.9 and also in three signal-to-noise-ratio, namely 6 < SNR < 8, 
8 < SNR < 10 and SNR > 10. The whole sky is divided in two different patch types (see, 
Fig. 1). Each type contains three different patches with different color codes provided in this 
table. 
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2.3 Error estimation 


In this paper, as we have mentioned earlier, we compare the statistical properties of the 
observed data in different direction of the sky without comparing with any theoretical 
model. To obtain the uncertainties/errors in the estimated PDF of the transmitted 
flux, we need to generate the covariance matrix associated with the flux PDF. Here 
too, we follow similar procedure as adopted in the BOSS analysis, in order to obtain 
the covariance matrix. We follow bootstrap resampling of the data in each bin. We 
have a total of 36 bins in both type of patch selection: 3 SNR bins, 3 redshift bins 
and 4 bins corresponding to each of the 3 patches and a complete sample. In each 
bin, we gather all the pixels that pass through our data cut and provided in Table 1. 
We perform bootstrap resampling of chunks upto lOOH wavelengths obtained from 
each quasar spectrum and generate 1000 realizations of the data. From these samples 
covariance matrix of the flux PDF can be easily calculated. Using different realizations 
we did check the convergence of the diagonal terms of the covariance matrix and we 
could conclude that the choice of 1000 realizations has been conservative. 

We concentrate on the statistical properties of the PDF or the data. Hence we 
calculate different statistical moments such as mean, median, variance, skewness and 
kurtosis of the PDF. We report the moments of the distribution that we obtain from the 
original data as well as the upper and lower error bounds on each of the moments that 
we calculate from the bootstrap simulations. For distribution of each moment, we gen¬ 
erate the empirical cumulative distribution function (ECDF) and report the (34.15%) 
upper and lower error bounds, which we hereafter shall refer as Icr uncertainties in the 
distribution of the particular moment. 

Before discussing the isotropy test on patches in next section, in Table 2 we report 
the mean transmitted flux and la uncertainties of the flux PDF that we obtain from 
the complete sample in the three redshifts and in three SNR bins. The number of 
quasars and the number of pixels that are used in our analysis are provided in Table 1 
for the complete sample size. Note that as we go back in time (higher in redshift) 
the increase of neutral hydrogen is evident from the decrease of the transmitted flux 
in all SNR. In each SNR, the error on the flux are comparable and that indicates, 
though at high SNR the signal is better, the less number of quasars in high SNR keeps 
the uncertainties comparable to the uncertainties at lower SNR which contain larger 
number of quasars (for example, compare the good and the best SNR cases). 

3 Analysis, results and discussions 

In the previous section we mentioned that we trisect the survey area in two ways. 
One patch, that is in southern galactic hemisphere remains common to both the selec¬ 
tion. For the 5 patches, that are color coded (red, green, blue, magenta and black) we 
obtain the PDF of transmitted flux. However, we compare the statistical properties of 
the PDF within each selection type. 

To start with, we provide the flux PDF with the errors in Fig. 2 for patch selection 
1 (corresponding to the selection in the left of Fig. 1) and in Fig. 3 for patch selection 
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Redshift range( 2 ;) 

SNR 

F±AF 

2.15-2.45 (z = 2.3) 

6-8 

8-10 

> 10 

0.826^°;^^^ 

0.822l°i[i« 

0.010_Q 4g7 

2.45 -2.75 {z = 2.6) 

6-8 

8-10 

> 10 

u. mz_Q 39 

U. <0O_Q 427 

U. 1 OU_Q 454 

2.75 -3.05 {z = 2.9) 

6-8 

8-10 

> 10 

O.Q87+^oill 

0 .686l°il^ 


Table 2. The flux statistics of the Lyman-a forest transmitted flux in different redshift 
bins and different signal-to-noise bins, provided for the samples across the complete sky. The 
mean flux and the Icr error bars on the PDF on both sides of the mean are provided. 

2 (the right selection in the Fig. 1). The colors of the PDF’s correspond to different 
patches (as shown in the Fig. 1). The PDF’s are calculated in the flux range [-0.2- 
1.5] in 34 bins. Theoretically the transmitted flux should be within 0 and 1 (for the 
complete and no absorption respectively) but due to the noise in observations and the 
bias in the continuum estimations, the transmitted flux might be obtained outside the 
theoretical boundary. The PDF’s are normalized such that the total area under each 
PDF is 1. In both the hgures, we provide the PDF for different redshifts and different 
SNR. As we mentioned in the previous section, due to less number of quasars in the 
higher SNR bins, the errors on the flux PDFs are comparable to that of the lower SNR 
bins. In the plots of the PDF, the difference between the PDF’s in good, better and best 
SNR bins are evident. For the higher SNR bins the flux PDF are more sharp around 
the peak. It is interesting to note that at the same SNR and at the same redshift, the 
PDF of transmitted flux from different patches are very similar. In some of the bins 
we hnd the flux from one patch is different from other ones. To quantify the difference, 
we next look for the different moments of the PDFs. 

In this paper, we restrict ourselves to calculating up-to the fourth moment, i.e. 
till kurtosis of the PDF. We calculate the mean [F), the median (F’ 1 / 2 ), the variance 
(cr^), the skewness (s) and the kurtosis {k) of the flux PDFs in each bin. Since we 
have noticed in Fig. 2 and 3 that the PDFs are very different in different redshifts, 
we report the residual of such quantities obtained in a patch w.r.t. the complete 
sample. In Fig. 4 we plot the 5 statistical quantities in residual space for the patch 
selection 1 and in Fig. 5 we plot the same quantities for patch selection 2. As we are 
plotting in residual space, we hnd the moments are distributed about the zero. We 
should mention that here the errors represent the uncertainties of the corresponding 
statistical moments obtained from the ECDFs of the moments. The plotted errors 
represent la uncertainties. 

Note that almost all the statistical moments from different patches are comparable 
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z = 2.3;6<S/N<8 


z = 2.3;8<S/N<10 


z = 2.3; S/N>10 



F 



F 

z = 2.9:6<S/N<8 



F 



Figure 2. Comparison of PDF of the Lyman-a transmitted flux in different patches for patch 
selection 1 (left of Fig. 1). The color codes represent the PDF’s from the corresponding patches. 
Along the rows we plot PDF’s for different redshifts and along the columns we plot for different SNR. 
The error on the PDF is estimated through bootstrap resampling over lOOA data chunks. 

to each other within their associated uncertainties. We would like to stress that we 
hud this consistency throughout redshift 2 — 3 and at all SNRs. Hence our results are 
consistent with the isotropic distribution of neutral hydrogen and also consistent with 
isotropic absorption of photons by the hydrogen clouds in the IGM at different redshifts. 

A closer examination shows that almost none of the patches contain any moment 
which is systematically higher/lower than the total sample during the time evolution. 
Hence the properties in each patch is preserved during the matter dominated expansion 
of the Universe. 

We should mention that we have obtained some deviations from isotropy in our 
analysis. We hud the residual moments for the blue patch and the red patch (for 
patch selection 1) do not agree at 2a level in cases. For example, we refer to the 
median mismatch (at z = 2.3, 6 < SNR < 8), the variance mismatches (at 2 ; = 2.9, 
6 < SNR < 8 and a.t z = 2.6, SNR > 10) and skewness mismatches (at 2 ; = 2.9 at 
8 < SNR < 10). Similarly, for patch selection 2 we hnd such deviations at 2a level 
within the red and the black patches. Since, the number of such deviations are only a 
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z = 2.3;6<S/N<8 



F 



F 


z = 2.9;6<S/N<8 



F 


z = 2.3;8<S/N<10 


z = 2.3; S/N>10 



Figure 3. Comparison of PDF of the Lyman-a transmitted flux in different patches for patch 
selection 2 (right of Fig. 1). The color codes represent the PDF’s from the corresponding patches. 
Along the rows we plot PDF’s for different redshifts and along the columns we plot for different SNR. 
The error on the PDF is estimated through bootstrap resampling over lOOA data chunks. 


few and all the moments agree within 3 ct, we do not report any statistically signihcant 
deviation. Moreover, we shonld note that dne to less nnmber of qnasars, the red patch 
at the sonthern hemisphere contains the largest dispersion in the distribntion of the 
statistical moments which in tnrn can lead to snch flnctnations. Hence, in order to rnle 
ont or conhrm any deviations between the bine and the red patches (or between black 
and red patches) we need more detection of Lyman-a forest in the sonthern galactic 
hemisphere. Fntnre surveys will be able to provide signihcantly more quasar spectra 
with higher sky coverage, which are essential to extend this initiative beyond the flux 
PDF level. 
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Figure 4. First few statistical moments of the Lyman-a transmitted flux PDF for different parts 
of the sky. For relative comparison we plot the residual moments from the total sample. The error 
bars correspond to the Icr bounds on that particular moment. This figure represents the comparison 
between the properties of the patches for patch selection 1 (left plot of Fig. 1). 
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Figure 5. First few statistical moments of the Lyman-a transmitted flux PDF for different parts 
of the sky. For relative comparison we plot the residual moments from the total sample. The error 
bars correspond to the Icr bounds on that particular moment.This figure represents the comparison 
between the properties of the patches for patch selection 2 (right plot of Fig. 1). 
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4 Conclusions 


As a first approach to test the isotropy in the matter dominated epoch we have used 
the Lyman-a forest data in this paper. In order to remain independent from theoretical 
assumptions of the IGM, we have used only the observational data and compared the 
statistical properties of the PDF of the observed flux in different directions of the sky. 
Detection of large number high redshift quasar spectra by BOSS has enabled us to 
perform this test. Though we report the distribution of neutral hydrogen is consistent 
with isotropic Universe during 2 ; ~ 2—3, making any general claim about the isotropy of 
the Universe requires much more sky coverage. As we have mentioned before, our test 
of anisotropy is partial due to only 3275 deg^ sky coverage. With a larger sky coverage 
we can address the precise direction and signihcance of the possible anisotropy, if 
present. With upcoming observations like e-BOSS [38] and DESI [39] we expect to 
detect higher quality of Lyman-a forest data, signihcantly higher number of quasar 
spectra with larger sky coverage. Hence, with the upcoming data the assumption of 
isotropy can be falsihed with higher precision. Cross-correlating the data from different 
surveys will be also important to rule out any systematic effect. 

Any presence of anisotropy in the Lyman-a forest is interesting and it points to¬ 
wards the distribution of neutral hydrogen, temperature-density relation and few other 
properties in the IGM. A straightforward extension of this topic would be to model the 
IGM using some semi-analytical modeling or simulations. With the modeling we can 
address if the signihcance of such anisotropy being statistical or physical. Moreover, if 
any physical anisotropy is found, we need to examine the change in the properties of 
the IGM it refers to and search for the probable cause. 

We would like to conclude by mentioning that the major hnding of our analy¬ 
sis is consistency of the data with the isotropic Universe in the final stage of matter 
dominated epoch (z ^ 2 — 3), and unfortunately we found no preferred direction in the 
Lyman-a forest to guide the travelers. Though our result is a partial test, if similar 
tests on larger survey area also conhrms the isotropy, that might hint towards the 
possibility that any late time anisotropy is probably caused by bulk how and not an 
intrinsic anisotropy in the Universe. 
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